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Abstract 

We describe an alternative approach to the prediction of W and Z transverse 
momentum distributions based on an extended version of the DDT formula. The 
resummation of large logarithms, mandatory at small q T , is performed in q T -space, 
rather than in the impact parameter b. The leading, next-to-leading and next- 
to-next-to-leading towers of logarithms are identical in the 5-space and g T -space 
approaches. We argue that these terms are sufficient for W and Z production 
in the region in which perturbation theory can be trusted. Direct resummation 
in g T -space provides a unified description of vector boson transverse momentum 
distributions valid at both large and small q T . 


1 Introduction 


We re-examine the transverse momentum distributions of vector bosons, in view of the 
large data samples expected at the Tevatron. In pp collisions at y/S = 1.8 TeV we 
expect about 10 5 W bosons and 10 4 Z bosons, observed through their leptonic decays, 
per 100 pb _1 of accumulated data. These events will be invaluable for QCD studies, as 
well as for precision measurements of the W mass. In order to exploit these data samples 
fully the experimenters will require detailed information about the expected rapidity and 
transverse momentum distributions of the vector bosons and of their decay products. 

In QCD a vector boson acquires transverse momentum q T by recoiling against one 
or more emitted partons §0- Order by order in perturbation theory we encounter 
logarithms, In Q 1 2 /q 2 , where Q is the mass of the lepton anti-lepton pair resulting from the 
vector boson decay. These logarithms must be resummed to give an accurate prediction 
in the low q T region. The original approach to the summation of logarithms at small q T 
was provided by Dokshitzer, Dykanov and Troyan (DDT) [j3| who derived an expression 
(reproduced here for the case of massive photon production), 

y = §E«^{AM(^.9 T )A7BK.? T )exp[T DDT ( ?T ,Q)] + fa » «)} . (1) 


where T is a leading log Sudakov form-factor,[] 
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The present state of the art in the theoretical description of vector boson production 
is based on the 6-space formalism where b is the impact parameter which is Fourier 
conjugate to the vector boson transverse momentum. The 6-space formalism, which allows 
the implementation of transverse momentum conservation for the emitted gluons,f] has 

1 Other notation will be defined in the body of the paper. 

2 See, for example, Ref. (J]. 
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the remarkable consequence that the cross section at q T = 0 is calculable for very large 

Q [0 §-0 

Nevertheless, in practice the 6-space formalism has certain disadvantages. Since the 
cross section is given as a Fourier integral in 6 which extends from 0 to oo, one cannot 
make theoretical predictions for any q T without having a prescription for dealing with 
the non-perturbative region of large 6. This problem can be solved by introducing an 
additional non-perturbative form factor (to be determined from experiment), but that 
also leads to unphysical behaviour of the cross section at large q T , where one should 
recover the ordinary perturbation theory result. These points will be further discussed 
later in the text. 

Clearly, if one could perform the Fourier integral in 6 analytically and thus obtain 
an expression for the cross section in g T -space, the above problems would be solved. 
A model for the non-perturbative region would have to be introduced only at the very 
lowest values of q T , and one would have a unified description of vector boson transverse 
momentum distributions valid at both small and large q T . 

In this paper we present an approach to resummation in g x -space, which is based 
on an extended version of the DDT formula. The 6-space formalism ||]-[|T|] resums the 
contributions to the cross section from the following towers of logarithms (L = In Q 2 /q^)- 


L : 

Qt 

NL : 

q-i 

NNL : 

Qt 

NNNL : 

Qt 


Our extended DDT expression agrees with the 6-space results for all but the NNNL series. 

3 Very massive vector bosons are produced at q T = 0 in association with semi-hard gluons which have 
zero net transverse momentum. Unfortunately, in W and Z production Q is too low and this asymptotic 
regime does not apply. 


3 



However, for vector bosons with masses less than Mz, we find that the NNNL series is 
numerically unimportant for q T > 3 GeV. Furthermore, a specific choice of coefficients in 
the g T -space Sudakov form factor allows us to absorb the first term in the NNNL tower 
of logarithms, and to obtain exact agreement with resummation in 6-space to 0{a 2 s ). 

Based on these results, we argue that the g x -space approach preserves almost all 
the reliable features of the 6-space formalism, [] and that it also has certain practical 
advantages: 


We avoid numerical pathologies in the matching, caused by combining results from 
6-space and g T -space. Although the matching is formally included in the 6-space 
method [10]. [T^. [T3J, the cross section is not correctly calculated for q T > Q/2. 
The cross section in this region is the result of a delicate cancellation between the 
resummed and finite pieces. The slightly different treatment of the two terms is 
sufficient to upset the cancellation. In contrast, the matching works well in q T - 
space, leading to a unified description of the q T and y distributions valid for all 
q T . 


• We need to introduce a model only at the very lowest values of q T . 


• We have the practical advantage that we avoid both the numerical Fourier transform 
and multiple evaluations of the structure functions at each value of q T . 

A complete explanation of these points will be found later in the paper. 

ft is important to emphasize here we are not challenging the theoretical importance of 
the 6-space formalism, which leads to interesting results, particularly about the production 
of very massive bosons. Nevertheless, it is our opinion that in practice the extended DDT 
approach is sufficient for the theoretical description of the W and Z production. 

The rest of the paper is organized as follows: in Section ^ we review the 6-space re- 
summation. In Section |3] we derive an extended version of the DDT expression, which 

4 The subleading terms in 6-space have a profound effect at q T = 0. However, for W and Z production 
this region is dominated by non-perturbative effects. 
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forms the basis of our approach. Section ff] contains comparison of the perturbative Su- 
dakov form factors in the g x -space and 6-space formalisms, and shows that in the region 
where the latter is reliable it is essentially identical with the former. We also present 
a prescription for dealing with the non-perturbative region of low q T , and compare our 
results with typical 6-space calculations. Our conclusions are given in Section || while 
Appendix |A| contains the saddle point evaluation of the 6-space expression for the cross 
section at q T = 0. 


2 Resummation formalism in 6-space 


The general expression for the resummed differential cross section for vector boson pro¬ 
duction in hadronic collisions may be written in the form 


da{AB -»■ V{-^ IV)X) 
dq% dQ 2 dy d cos 6 d<fi 


Q 2 


x 


2 8 NttS ( Q 2 - M 2 ) 2 + M 2 T 2 V 
Y r (q 2 T ,Q 2 , y ,9) + Y f (ql,Q 2 , y ,9,<t )) 


(4) 


In the above, N = 3 is the number of colors, \/S is the total hadron-hadron center- 
of-mass energy, while 9 and (j) refer to the lepton polar and azimuthal angles in the 


Collins-Soper (CS) frame []n|. The mass and width of the vector boson are denoted by 
My and Ty. The functions Y r and Yf stand for the resummed and finite parts of the cross 
section, respectively. As the details of the finite part are not important for the subsequent 
discussion, we review here only the resummed part, and refer the reader to Ref. |T3) for 
the complete description of 0(as) finite part. 

The resummed part of the cross section is given as the Fourier integral over the impact 
parameter 6, 5 

1 


Yr 


POO 

(qlQ\y,») = 0(Q 2 -qD— dbbMq T b)Y.’FS, F {QA x A ,x B ) 

Z7T JO 77 

a,b 


3 The prime on the sum in Eq. ([|) indicates that gluons are excluded from the summation. 
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bo, r, / b 0 . 


x W ab (Q,h,9)K /A (x A ,^)fl /B (x B ,^) , 


(5) 


bj 0/ ° K K 

where the variables x A and x B are given in terms of the lepton pair mass Q and rapidity 
V as 

n n 

( 6 ) 


Xa = exp ^ ’ Xb = exp 


The modified parton structure functions in Eq. ( 0 ), /', are related to the MS structure 
functions, /, by a convolution 
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where (a, b ^ g) 


, 7T~ 


C ab (z,n ) = <5 afe < ^(1 - z) + a s (/i)C F \l - z + (— - 4)5(1 - z) 


C a g(z,/i) = a s {ii)T R \2z(l - z) 
Here we have introduced 

&s{y) = 


®s(d) 

2vr 


(7) 

( 8 ) 

(9) 

( 10 ) 


while Cf = 4/3 and T R — 1/2 are the usual colour factors. 

The function W can be expressed in terms of the Sudakov form factor S(b, Q) and is 
given by 

WM, b, 6) = H f “\0) exp [S(b, (3)] , (11) 

where which includes the angular dependence of the lowest order cross section and 
coupling factors, is defined in Appendix A of Ref. 0. The Sudakov form factor itself is 
given as [6] 


S(b,Q) = - 


rQ 2 dfi 2 


'b 2 0 /b 2 /A 


111 V'iW'.n + B(asM) 


( 12 ) 


with bo = 2exp(—ye) ~ 1.1229. The coefficients A and B are perturbation series in as, 

OO OO 

A(a s ) = E i B(a s ) = E • (13) 


i= 1 


i— 1 


6 














The first two coefficients in the expansion of A and B are known ||T6l IT? 


AW 

= 2 C F , 
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5 (1) 

= —3 Cf , 
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(14) 


One of the main advantages of the 6-space resummation formalism is that the simple 
form for S(b, Q) as given in Eq. (]T2|), remains valid to all orders in perturbation theory [6j]. 
In addition, as mentioned above, for very large values of the vector boson mass the 6-space 
formulae make definite predictions for the q T = 0 behaviour of the cross section || j(|. 

Unfortunately, the practical implementation of the 6-space formulae presents some 
difficulties. The 6-space integral in the Bessel transform in Eq. (^) extends from 0 to oo, 
which means that one has to find a way to deal with the non-pertnrbative region where 6 
is large. That problem is usually circumvented by evaluating W and the parton structure 
functions at 


Vi + (&/<> Um) 2 ’ 


(15) 


which never exceeds the cut-off value 6 i im , and also by introducing an additional function 
F np , which represents the non-pertnrbative (large 6) part of the Sudakov form factor, 
to be determined from experiment 0. This is usually done by assuming a particular 
functional form for F NP which involves several parameters that can be adjusted in order 
to give the best possible description of experimental data. The specific choice of the 
functional form for F NP is a matter of debate i, MM, but we will not discuss it further 
here. The point which we would like to emphasize here is that without introducing 6* 
and F NP one would not be able to make theoretical predictions for any value of q T , even 
in the large q T region where perturbation theory is expected to work well. 
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Another problem which occurs in the 6-space resummation formalism is the transition 
between the low and the high q T regions. At large q T the resummed part is well represented 
by the first few terms in its perturbative expansion. When the resummed part Y r is 
combined with Yf one formally recovers the perturbation theory result. However, the 
cancellation at large q T is quite delicate and is compromised by the non-perturbative 
function which acts only on Y r . We illustrate the problem in Figure [T].[| which compares 
the 0(as) perturbation theory result for da/dq T in W + + W~ production at the Fermilab 
Tevatron, to the theoretical prediction obtained from the 6-space resummation (Eqs. (|4|) 
and ©)•□ 

Even though by carefully matching the low and high q T regions one can reduce theo¬ 
retical errors and produce smoother transverse momentum distributions, matching is still 
bound to fail eventually, and one is forced to switch to the pure perturbative result at 
some q T [|10]. This procedure inevitably leads to discontinuous q T distributions, which are 
clearly unphysical. 

If one could find a g T -space expression for Y r , both of the above problems would be 
solved: just as for the conventional perturbation theory, theoretical predictions could be 
made without any smearing or additional functions, at least for values of q T not too close 
to zero. Also, since Y r and Yf would both be calculated in g T -space, the cancellation 
between the resummed part and subtractions from the finite part would be explicit, and 
matching of Y r + Yf onto the perturbative result at large q T would be manifest. With this 
motivation we consider the derivation of g T -space equivalent of Eq. © in the following 
section. 


6 Note that throughout the paper we use the MRSR1 structure functions, with as(Mz) = 0.113 jl8|| . 

7 Following Eqs. (22,23) of Ref. ]li|, instead of & 0 /6* we actually used the exact first order result for the 
scale at which parton distribution functions are evaluated. This prescription preserves the total integral 
and reduces to &o/6* for large b. Furthermore, it improves the large q T matching between Y r + Yf and 
the perturbation theory. 



3 Resummation in g T -space: extended DDT formula 


For the sake of simplicity we discuss only the resummed part of the non-singlet (NS) cross 
section for the process AB —> 7 *X . The extension to the general process AB —> V(—> 
ll')X is straightforward. I 11 this case Eqs. (£D, (|5|) and (pT|) can be rewritten in the form 

da ~ 


°o 
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dq\ dQ 2 




f dx A dx B 5(x a x b - %-) 
Jo o 
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2 Jo 


dbbJ 0 (q T b ) exp [5(6, Q)] f q/A (x A , j-) f' q/B (x B , , (16) 


where a 0 = 4yra 2 /(95) and f' q/A = f' q/A - f q/A , f' q/B = f' q/B - f' q/B are the higher order NS 
structure functions. Note that we have removed the non-perturbative function F NP and 
variable 6 * from Eq. (|5|), so that the above expression represents the pure perturbative 
result. 

From Eq. (|T 6 | ) one can easily obtain the iV-th moment of the cross section with respect 
to r = x A x B = Q 2 /S, 


m) 



Q 2 da 

a 0 dq B dQ 2 


1 /»00 h h 

£ i db b Jo{qTb) exp Q)] tU(x j ) 4bW f) ■ 


(17) 


The IV-th moment of the NS higher order structure function satisfies the GLAP equation,Q 

= ) ( 18 ) 

with the solution 


60 . 


fq/H( N , j) = exp 

Using Eqs. (p~7| , [T9|) we may write 


rQ 2 

/(6o/b)2 /i 


^-7w0s(/^)) 


fq/A( N ,Q) ■ 


E(N) = G(N,Q)l r dbbJ 0 (q T b) exp[U N (b, Q)] , 
Z Jo 


(19) 


( 20 ) 


8 The anomalous dimension 7' differs in a calculable way from the MS anomalous dimension. 
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where G(N, Q ) denotes the parton flux, 


G(N,Q) = E4 f'„A(N,Q)fi /B (N,Q) , 

9 

and the exponent U is given as 


( 21 ) 


+v + Q) — — [ 
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bg/fe 2 fj, 
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Here iD 2 = —etc. Inserting Eq. (fZ2| ) in Eq. (^0[) we obtain 


S(iV) = G(N,Q) 


1 r 00 


db b J 0 (q T b) exp 


00 n+1 


E E aSW) i"” 

.n=l m=0 


' Q 2 b 2 ' 


bl 


n D r ri 


This expression may be integrated by parts using the relationship 




= xJn(x) 


( 22 ) 


(23) 


(24) 


Because of the rapid damping of the Sudakov factor as b —> 00 we may ignore the boundary 
terms and obtain 

1 __ f°° . d 


) /•OO W 

s(/v) = ~2A G(n ' q) L Wl(l) 5 txI 


00 n+1 / n2„2 ' 
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/•OO r] 
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'QV' 
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(25) 


Eq. (g5D already has the structure of the DDT formula. In fact, setting ln+feo = 0 in the 
integrand we recover exactly the DDT formula, i.e. the exponent has exactly the form of 


Eq. (22) with b^/b replaced by q T . Because of that we write E(7V) in the form 


d f r°° 

E(iV) = —+ (j(iV, Q) / dxJH+exp 
ttO” ' Jo 
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where the remainder R is defined as 


R(q T ) = G(N,Q ) J dx J\ (x) < exp 


OO n+1 

E E « 3 (< 3 ) In" 

_n=l m =0 


2 ^ 2 ' 


Q 2 x 
, q.T h l 


rD r 


— exp 


oo n+1 / 

E E “S(« i» m ' 

.n =1 m=0 


ql 


,D r 


( 27 ) 


Usingf] 


f°° dx R(x) {l, In f, In 2 f, In 3 f,... } = {l, 0, 0, - + (3),... } , (28) 

Jo 1 6 0 6 0 6 0 > y 2 > 

we can evaluate R(q T ) as a power series in as- We find that the remainder contributes to 

the NNNL tower of terms, three logarithms down from the leading terms (L = In Q 2 /q 2 ), 


with 


R(Qt) = ~G(N, Q)\ C(3) E rj{iD 2 ) j a j s {Q)L 2j ~ z + 0[a 3 s L 23 
l j=2 


r ) fo 40 28 , 11 13 , 

+ 2 , r 3 , r 4 , r 5 , r 6 , r 7 ,... j = |8, —, —, 4, — j . 


—4 \ 


(29) 


(30) 


3 3 9 45 

Starting from the 6-space expression we have demonstrated an extended DDT formula, 
da - n2 


dq 2 T dQ 2 


<Jq_ 

Q 2 

d 


E e 2 f dx A dx B <5+ A x B - %) 

a JO b 


x ^ j f q / A ( x A, ? T ) f q / B ( x B, g T ) exp [T(g T , Q)] + 0(a 3 s L 2j 3 ) j , (31) 

which holds if we drop NNNL terms. In the above expression the g T -space Sudakov form 
factor is given by 


T(q T ,Q) = - 


rQ 2 djx 2 


A{a s {n)) In ^ + B(a s (n)) 
d 


(32) 


where the g T -space coefficients A and I? are defined in a similar way as their 6-space 
counterparts, i.e. 

OO OO 

+«s) = E c+i (i) , B(a s ) = E ■ (33) 

2—1 2—1 

9 See, for example, Ref. |J]. 
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The first two coefficients in A and B would be exactly the same as corresponding 6 -space 
coefficients if we drop NNNL terms. However, by making the particular choice of 


A® = aw , 

^( 2 ) = A (2) 

bW = H (1) , 

BW = H (2) + 2(H (1) ) 2 C(3) , (34) 


we absorb the first term in the NNNL tower of logarithms. In this way Eq. ([54|) imposes 
exact agreement between the 6 -space and g T -space formalisms at order a 2 s . 

As we pointed out at the beginning of this section, the extension of the NS cross 
section for AB —> 7 *X to the general case of AB —» V(—> ll')X which includes the decay 
presents no difficulties, so that our g T -space equivalent of Eq. (||) is given in the extended 
DDT form as 


Y r (ql,Q 2 ,y,9 ) = Q(Q 2 - g 2 ) 


x 


Y,' H ab{ e )j^ [f'a/A^ A ,qT) K/b^b^t) exp[T(g T ,Q)]J , (35) 

a,b 


71 

d 


dql 


with T given in Eq. (^) in terms of coefficients of Eqs. ([k|,[||). The above equation is 
the central result of this paper. It is still ill-defined in the small q T region, which reflects 
the fact that the problem is not entirely determined by perturbation theory and requires 
non-perturbative input. We will discuss our model for the non-perturbative region later 
in the following section. 


4 Results 

4.1 Form factors 

Before presenting our results for W and Z production we compare the form factors cal¬ 
culated using the 6 -space and g T -space formulae, for values of Q which are presently of 
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interest. In practice this means Q < M Z - The comparison of the form factors will allow 
us to make an estimate of the practical numerical importance of transverse momentum 
conservation, i.e. of the subleading terms which are not present in the g T -space formalism. 
To simplify the comparison we will consider the effects of the Sudakov form factor alone. 
We will therefore ignore the influence of modified parton distribution functions on the q T 
dependence. For the purpose of illustration we take Q = M z and as(M z ) = 0.113. 

We define the 6-space form factor as 

Q 2 


F {b \q T ) = J d 2 b exp {ib-q T ) exp [S(b*,Q)] F NP (Q,b,x A ,x B ) 

rp /.oo 

= — dbbJ 0 (bq T ) exp [S(K,Q)} F NP (Q,b,x A ,x B ) . 

Z Jo 


(36) 


Note that F NP and 6* have to be introduced in the above expression as a prescription 
for dealing with the non-perturbative region of large b. A specific choice of the non- 
perturbative function should make a difference only in the region of low q T . In order to 
show that, in Figure [2| we present form factors evaluated with F NP taken from Ref. JTT 

(ERV, g = 3.0 GeV 2 


(LY), and with an effective gaussian as used in Ref. |13[ (ERV, g = 3.0 GeV z ). For LY 
form factor we take x A = x B = M z /\fS for \/S = 1.8 TeV. As expected, at large q T 
the form factors resulting from the two choices of F NP agree. For small q T we find that 
results for F ( F (q T ) tend to a different finite intercept controlled by the non-perturbative 
function. 

The above 6-space expression for the form factor should be compared to its g T -space 
counterpart, 


d 


F {qT) (q T ) = Q 2 —irexp[T(q T ,Q)} , 


dq% 


and also to the d(a|) perturbation theory result, 

r\1 2 2n—1 r \2 


Qt n =1 m =0 

with n C m given in terms of the g T -space coefficients as 

iC 1 ! = A (1) , 


'■It 


(37) 


( 38 ) 
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iC 0 = B « , 

2C3 = -^(i (1) ) 2 , 

2 C 2 = i«(/3 0 -|sW), 

2Cf = AM + bW(Po-bW) , 

2 C 0 = 5 (2) . (39) 

As one can see from Figure |3|, in the region where one can trust perturbation theory 
(q T > 3 GeV), our g T -space result of Eq. (|37|) agrees well with the 6-space form factor 
(obtained with ERV non-perturbative function). Further, it is clear that resummation is 
needed in the region where F ^ and F^ differ significantly from the perturbative result. 

It is also interesting to investigate the size of the NNNL effects. In Figure |]we show 
the g T -space form factor F^ T '(g T ) calculated using coefficients given in Eq. (|34|), and also 
the one calculated with B ^ replaced by B^ 2 \ As one can see, the change is never more 
than a few percent for q T > 3 GeV. 

We therefore conclude that the 6-space and g T -space formula are substantially identical, 
despite the neglect of NNNL terms in the latter. The differences between them are smaller 
than the differences introduced in the 6-space formalism by the use of different non- 
perturbative functions. The above conclusion holds for the particular case of the vector 
boson production with Q < M z . 


4.2 Extension to the non-perturbative region 


As we have already pointed out, there are two main advantages of the g T -space approach 
over the 6-space formalism: first, outside of the non-perturbative region one can make 
theoretical predictions based on perturbation theory alone with soft gluon resummation 
effects included. In Figures |] and |(| we show predictions of Eq. (|35|) for W + + W~ and 
Z production at Fcrmilab Tevatron. It is clear that these predictions are quite close to 
typical 6-space results. Second, matching of the resummation formalism onto pure pertur- 
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bation theory for large q T is explicit, and hence there is no need for somewhat unnatural 
switching from one type of theoretical description to another. Since our calculation con¬ 
tains the 0(as) finite part and the 0(a |) Sudakov form factor, there still may be some 
residual unmatched higher order effects present in da/dq T in the large q T region, where the 
cancellation of the resummed part and subtractions from the finite part is quite delicate. 
However, these effects are expected to be small, and should be even less important after 
the inclusion of the second order calculation of Yf. The g T -space matching is illustrated 
in Figure |7] for W + + W~ production at Tevatron, and should be compared to the 6 -space 
result shown in Figure [lj. Note that less than 2% of the total cross section lies above 
q T = 50 GeV, so the overall importance of the portion of the cross section shown in Fig. [7| 
is quite small. 

Up to now we have discussed only the g T -space predictions in the perturbative region, 
i.e. for q T > 2 — 3 GeV. Still, in order to compare theoretical predictions to experiment 
one has to find a way of dealing with the non-perturbative region (q T —> 0), where Eq. (j35|) 
is ill-defined. The form of Y r suggests that we make the following replacement in Eq. (|H5|) : 


fa/A&A, Qt) fb/ B {x B, q T ) exp [T(q T , Q)] — > 

fa/ A (Ta.) Qt*) f' b / B { x B; (It*) exp [T(g T *, Q)} F NP (q T ) . (40) 

Here, g T * is the effective transverse momentum and F NP is the g T -space non-perturbative 
part of the form factor. Since the above replacement should affect only the region of small 
q T , we define g T * as 

r ^ 1 (41) 


Qt* (It T ^TUm exp 


Qt 

^Tlim 


which never goes below the limiting value g T iim, and also approaches q T as q T becomes 


dql 


F np we require that 




F NP (q T ) - 

+ 0 

(for 

q T 

-0) , 

F NP (q T ) - 

-»• 1 

(for 

q T - 

+ Q) 

F NP (q t) - 

■+ const. 

(for 

q T - 

-o). 


(42) 
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The first two properties ensure that the integral of Y r over q 2 gives the result 


rQ 2 




( 43 ) 


a, b 


which is required to reproduce the exact O(as) total cross section after integration over 
q T , as explained in Ref. |T3j. The third condition is motivated by the analytic 6-space 
results for da/dq 2 in the limit where q T —► 0 || (see Appendix [A|). 

A simple choice for F NP which satisfies all of the above requirements is 


F NP (q T ) = l~exp[-aql] 


(44) 


At q T = 0 this yields 


da 

dq 2 

d / da 


oc a iP) [fa/Ai X Al 'Alim) fb/B ("^B ) Alim) eX P [^(^Tlim) Q)] 


a,b 


dq 2 ^ _ ® 2 E [/«/i( lA > fe-) /vb( X b > ? Tli J exp [T(q Tlim , Q)]] .(45) 

Therefore, a and g Tlim control the intercept and the first derivative of da/dq 2 at q T = 0. 
The effects of changing these non-perturbative parameters are illustrated in Figures and 
H for W + + W~ production at Fermilab Tevatron. In Figure |8| we compare typical 6-space 
results for the da/dq 2 distribution, to the g T -space predictions with several different values 
of a (g-riim was fixed to 4.0 GeV). In Figure |9| we plot our da/dq T results obtained with 
a fixed to 0.10 GeV~ 2 , and for several different choices of g Tlim . These results show how 
varying g Tlim modifies the width and shifts the peak of the da/dq T distribution. 

From Figures Rl and 0 it is also clear that a and g Tlim affect only the low q T region, 


while for q T > 10 GeV we again obtain the extended DDT result of Eq. (|35|). Because 
of that, determination of these parameters from the experimental data should not be too 
difficult. 
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4.3 Overall smearing 


Introduction of a and q Thm allowed us to extend the validity of Eq. (|35|) beyond the 
perturbative region in q T . However, this may not be enough for a good description of 
experimental data, and one may need additional degrees of freedom for modelling the low 
q T region.^ This can be achieved by choosing more complicated functional forms for g T * 
and F NP than the ones we suggested in Eqs. (41, 44), or by imposing an overall smearing 
on the theoretical transverse momentum distributions. Here we briefly discuss the later 
possibility. 

Suppressing irrelevant variables, the smeared cross section is given in terms of 


YMl) = J d 2 k T f( |k T - q T |) Yi{kl) , (46) 

where Y % stands for either resummed or hnite part in g x -space, and / is the smearing 
function. For the sake of simplicity we take a gaussian, 

f(k T ) = -exp(-<)/4) , (47) 

7T 


with g being an additional non-perturbative parameter. The above choice is convenient 
since the azimuthal integration can be done analytically. This leads to the final expression 
for the smeared Y), 


Yi{ql) =g dkl exp - g(ql + kl) I 0 (2gq T k T )Y,(kl) 


(48) 


Note that both the resummed and the hnite part of the cross section can be smeared 
together, and therefore the smearing procedure should not affect matching onto the pure 
perturbative result at large q T . The effects of an overall smearing for low q T are illustrated 
in Figure [L(J for W + + W~ production at Tevatron. 

10 We remind the reader that some choices of the non-perturbative function in the 6-space formalism 
involve 4-6 different parameters. 
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5 Conclusions 


In this paper we have outlined an approach to the calculation of the transverse momentum 
distributions of W and Z bosons using an extension of the DDT formula which works 
directly in g T -space. Our formalism agrees with 6-space for all calculated logarithms except 
the NNNL series. This is a pragmatic approach which uses the available theoretical in 
an efficient way. For q T above about 3 GeV the cross section is essentially determined by 
perturbative QCD. In the region q T < 3 GeV the cross section is determined by a model, 
the form of which is motivated by the analytic results from the 6-space approach. Just 
as in the 6-space approach, the details of the model are to be fixed by comparison with 
experiment. The numerical program incorporating our results describes all kinematic 
regions. 

An obvious shortcoming of this paper is the failure to include the results of the order 
a 2 s calculations JT?|, [T9], ^Cj (generalized to include the decay of the vector boson 0.0) 
in the finite part of the cross section. In the g T -space formalism these should be relatively 
straightforward to include. After inclusion of these effects we will have a full description 
of vector boson production valid in all kinematic regions, with a minimum of model 
dependence. 
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A Analytic behaviour at q T = 0 

The result of Parisi and Petronzio 0 for the intercept at q T = 0 can be obtained by saddle 
point evaluation of Eq. (|5E|), 



( 49 ) 
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In writing the above equation we have assumed that the saddle point value of b is small 
so that 6* = b and F NP (b) = 1. Introducing the variable x = ln& 2 we have that 

F ( - b \0 ) = J efeexpf— h(x)] , (50) 

where 

h{x) = -[x + 5(exp(x/2), Q)\ . (51) 

The saddle point result for F^ b \ 0) is then given by 

F {b) (0) = -F-J r^ 77 \ exp[-/i(ar SP )] , (52) 

4 V h"{x SP ) 

where a: SP = In 6 2 p is defined by the condition 

h\xsp) = 0 . (53) 


On the assumption that the structure functions are slowly varying functions of the scale, 
the resummed part at q T = 0 becomes 


Y r (0,Q 2 ,y,6) 


& 2 p I 2n 
47T V -S"(b SP ,Q) 


^2'W ab (Q,b SP ,d)f^ /A (x A , 

a,b 


bo_ 

bsp 


)fb/B^ X 


B 1 



(54) 


where 


S"(b,Q) = 


d 2 S(b,Q) 


(55) 


d(ln & 2 ) 2 

By retaining only the leading term (Ai) in the Sudakov form factor we can obtain 
an approximate analytic solution. We assume that the running coupling satisfies the 
equation {(3q = (33 — 2n/)/6) 

2vr 1 


«s(/ i ) — 


(3q lnfi 2 /A 2 ’ 


(56) 


and set C = 2Cp/f3o. The saddle point of the integral is then given by 

j_ = A 

b s P b 0 'A' ' 


(57) 
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Using Eqs. (|54| , |57|) we obtain the final result for the resummed part of the cross section 
(L = lnQ 2 /A 2 ), 


Y r (0,Q 2 ,y : d) 


h 2 
°0 


\/2ttCL A 2 \v 


4ttA 2 (C + 1) V Q 2 


AYY,' H Smf«/A( x *’T L KB(x»’T L ) . ( 88 ) 


a, b 


with 

r/ = C In + 1 . (59) 

The Y) has a finite intercept at q T = 0 which shrinks with Q. For Uf — 3 (4) we have that 
ri = 0.586 (0.602). 
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W + +W production in pp at VS=1.8TeV 
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Figure 1: Comparison of the 6-space da/dq T distribution for W + + W~ production at 
a /S = 1.8 TeV with 0(as) perturbative calculation. The resummation results were ob¬ 
tained with pure gaussian (g = 3.0 GeV 2 .6ii m = 0.5 GeV -1 ) form of F NP . We assumed 
BR(W —> ev) = 0.111. 
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F(q T ) 



Figure 3: Form factors F lyb \ F^ and F^ p \ The 6-space results were obtained with an 
effective gaussian form of F NP (g = 3.0 GeV 2 , b n m = 0.5 GeV -1 ). 
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Figure 5: Comparison of various theoretical predictions for W + + W~ da/dq T with CDF 
data |^3[. The 6-space results were obtained with an effective gaussian form of F NP 
(g = 3.0 GeV 2 ,6 lim = 0.5 GeV -1 ). We assumed BR(W —*• ev) = 0.111. 
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d<r/dq T [pb/GeV] 



Figure 6: Comparison of various theoretical predictions for Z da/dq T with CDF data 
fhf|j . The 6-space results were obtained with an effective gaussian form of F NP (g = 
3.0 GeV 2 ,6i im = 0.5 GeV -1 ). We assumed BR(Z —»■ e + e“) = 0.033. 
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Figure 7: Comparison of the q T -space da/dq r distribution for W + + W production at 
a fS = 1.8 TeV with 0(as ) perturbative calculation. We assumed BR(W —» ev) = 0.111. 
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Figure 8: Various theoretical predictions for da/dq 2 in W + + W~ production at \/~S = 
1.8 TeV. The 6-space results were obtained with an effective gaussian form of F NP (g = 
3.0 GeV 2 , b i; m = 0.5 GeV -1 ). The g T -space predictions correspond to g Tlim = 4.0 GeV. 
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Figure 9: Various theoretical predictions for da/dq T in W + + W~ production at \/~S = 
1.8 TeV. The 6-space results were obtained with an effective gaussian form of F NP (g = 
3.0 GeV 2 , b i; m = 0.5 GeV -1 ). The g T -space predictions correspond to d — 0.10 GeV -2 . 
We assumed BR(W —> ev) = 0.111. 
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Figure 10: Effects of smearing in g T -spaee for W + + W production at \/~S = 1.8 TeV. 
We used a = 0.10 GeV -2 , g Tlim = 4.0 GeV, and BR(W —> ev) = 0.111. 







